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p i' Abstract 



The Newtonian solid-mechanical theory of non-compressional spheroidal and tor- 
sional nodeless elastic vibrations in the homogenous crust model of a quaking neutron 
r~| I star is developed and applied to the modal classification of the quasi-periodic oscil- 

lations (QPOs) of X-ray luminosity in the aftermath of giant flares in SGR 1806-20 
P I and SGR 1900-1-14. Particular attention is given to the low-frequency QPOs in the 

c/5 , data for SGR 1806-20 whose physical origin has been called into question. Our calcu- 

lations suggest that unspecified QPOs are due to nodeless dipole torsional and dipole 
^ ! spheroidal elastic shear vibrations. 
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1 Introduction 

This study was undertaken as a part of current extensive theoretical investigations of aster- 
oseismology of a neutron star[Tl [21 El IH O [6] (see also references therein) that have been 
boosted by recent discovery of quasi-periodic oscillations (QPOs) of X-ray luminosity in the 
decaying flares of two magnet ars [71 E], SGR 1806-20 and SGR 1900-1-14, with concomitant 
suggestion to interpret this variability as caused by quake induced differentially rotational, 
torsional, oscillations. Following this suggestion the focus of most theoretical works is on 
computing the frequency spectra of odd-parity torsional mode of shear vibrations and less 
attention is paid to the even parity spheroidal elastic mode. However, from the viewpoint of 
modern global seismology [91 ITO] . the spheroidal vibrational mode in a solid star and planet 
has the same physical significance as the toroidal one in the sense that these two funda- 
mental modes owe their existence to one and the same restoring force [Tl] IT2l 113] . In this 
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light there is a possibihty that, by not considering both these modes on an equal footing, we 
may miss discovering certain essential novelties which are consequences of solid mechanical 
laws governing seismic vibrations of superdense matter of neutron stars. Adhering to this 
attitude and continuing our current invest igat ions |13j, we derive here spectral equations for 
the frequency of both spheroidal and torsional elastic nodeless vibrations in the solid crust of 
quaking neutron star and examine what conclusions can be drawn regarding low-frequency 
QPOs whose physical nature still remain unclear. 

2 Frequency of nodeless spheroidal and torsional elas- 
tic shear vibrations in homogeneous crust 

In this paper we follow the line of argument of the standard two-component, core-crust, 
model of quaking neutron star[Tl] in which crustal metal-like material (composed of nuclei 
dispersed in the sea of relativistic electrons) is treated as a highly robust to compressional 
distortions elastic material continuum of a uniform density p characterized by constant value 
of shear modulus /i. In this model it is presumed that the quake-induced non-compressional 
seismic vibrations driven by bulk force of pure shear elastic deformations (which are not 
accompanied by fluctuations in density 5p = —pVkUk = 0) can be adequately modeled by 
equation of Newtonian, non-relativistic, solid mechanics 



From now on ^^(r, t) stands for the field of material displacements in the crust of the depth 
AR = R — Rc with R and Rc being radii of star and core, respectively. The linear relation 
between tensors of shear elastic stresses and shear deformations or strains Uik is the 
Hooke's law of elastic (reversal) shear deformations. 

In what follows we focus on poorly investigated regime of nodeless shear vibrations in 
which the fields of oscillating material displacement subject to the vector Laplace equation [13] 



This last equation can be thought of as the long wavelength limit of vector Helmholtz 
equation describing standing- wave regime of vibrations[51 |T3]. In these latter works it has 
been shown that the eigenfrequency problem of global nodeless vibrations of a solid star can 
be unambiguously solved with aid of the Rayleigh's energy method. In the present paper 
this method is extended to the case of spheroidal nodeless vibrations which are considered 
in one line with torsional ones. 

The stating point of the energy variational method is the integral equation of the energy 
balance 



pui = VfcO-ifc, aik = 2puik, Uik = - 



[ViUk + VfcMj], Ukk = ^kUk = 0. (1) 




(2) 
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which is obtained by scalar multiphcation of equation of sohd mechanics, ([T]), with Ui and 
integration over the volume of seismogenic layer. The field u(r, t) can be conveniently rep- 
resented in the following separable form 

u(r,t) = a(r)a(t) (4) 

with the field of instantaneous (time-independent) displacements a(r) obeying, as follows 
from ([2]), to equations 

V2a(r) = 0, V ■ a(r) = (5) 

and a{t) stands for the temporal amplitude of vibrations. Inserting (jll) in ([3]) we arrive at 
equation for a{t) having the form of standard equation of normal oscillations 

dE ^ ^ Ma^ Ka'^ .. 2 ^ 2 K 

_=0, — + — ^ „ + .V = 0, = (6) 

M = J pttittidV, K = 2 j fxaihttikdV aik = ^[ViUk + V kai\. (7) 

The solenoidal fields of instantaneous material displacements in two fundamental modes of 
nodeless vibrations - the spheroidal (normally abbreviated as oSi) and the toroidal (ab- 
breviated as 0^^), are determined by two fundamental (orthogonal and different in parity) 
solutions to the vector Laplace equation which are uniquely defined by the general solution 
to the scalar Laplace equation V^x(i') = 0- spherical coordinates with fixed polar axis, 
the solution of (I5l) corresponding to nodeless spheroidal vibrations, a^, is given by the even 
parity poloidal (polar) vector field and instantaneous displacements in the torsional mode, 
a^, are described by the odd parity toroidal (axial) vector field: 

a, = Vx Vx (rx), a^ = V x (rx), (8) 
X(r) = /,(r)P,(cos^), /,(r) = [Ae/ + (9) 

Henceforth P^(cos^) stands for the Legendre polynomial of multipole degree £ and and 
are the arbitrary constants to be eliminated from boundary conditions on the core-crust 
interface and on the star surface. Thus, in order to obtain the frequency spectra of both 
spheroidal and torsional nodeless vibrations one need to specify first these constants and 
then to compute integrals for integral parameters of vibrations, that is, the inertia M and 
the stiffness K. Some useful mathematical details of such calculations can be found inllSI. 



2.1 Spheroidal mode 

The poloidal field of nodeless instantaneous displacement a^ in s-mode is irrotational: V x 
a^ = 0. To specify Ai and we adopt on the core-crust interface, r = Rc, the condition 
of impenetrability of seismic perturbation in the core. On the star edge, r = R, we impose 
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the condition that the radial velocity of material displacements equals the rate of spheroidal 
distortions of the star surface [T5l 



r\r=Rc 



0, 



u 



r\r=R 



R{t), R{t) = R[l + a{t)Pi{cost 



(10) 



The solution of resultant algebraic equations leads to following values of arbitrary constants 



Ai 



+ 1)' " ' i?2m_i?2m- 

Tedious but simple calculation of integrals for inertia M and stiffness given by 
poloidal field yields 

A'kR^P 
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(11) 
with 

(12) 

(13) 
(14) 



The fractional frequency of nodeless spheroidal irrotational shear vibrations as a function of 
multipole degree I is given by 



00. 



2(2£+ 1) 
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(£2 - 1)(1 - A2^-l) + i{e + 2)A2^-l(l - A2^+3) 



(l-A2^+i) _ 
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(£ + 1) + £A2^+i 
t^s(^) = 27rz/(oS^)]. 



(15) 
(16) 



It worth noting that in the limit of zero-size radius of the core, A = (Rc/R) —>■ 0, when 
entire volume of the star sets in vibrations, we regain the early obtained spectral formula 
for global nodeless spheroidal nodeless shear vibrations [13] i^ioSi) = z/q [2{2i + 1){£ — l)]^^'^ 
showing that the lowest overtone of the global nodeless spheroidal oscillations in the entire 
volume of the star is of quadrupole degree, £ = 2 (see alsojlT]). In the meantime, the lowest 
overtone of spheroidal vibrations trapped in the crust is of the dipole degree, £ = 1. This 
suggests that the dipole overtone can be considered as a signature of spheroidal vibrations 
locked in the crust. It is remarkable that dipole vibration can be thought of as, so called, 
Goldstone's soft mode whose most conspicuous feature is that the frequency as a function 
of intrinsic parameter A of oscillating system uJs{£ = 1, A) ^ 0, when A ^ 0. In the model 
under consideration this parameter is given by A = (Rc/R). The limit A = belongs to 
translation displacement of the center-of-mass of the star, not a vibration; this is clearly 
seen from the equation for energy (Hamiltonian) of harmonic oscillations (EI). 

In the upper panel of Fig.l we plot the fractional frequency Us{£)/uq as a function of 
h = AR/R. In the next section we show how the input parameters of obtained spectral 
equation (fTSll . namely, the natural unit of frequency z/q and the depth h of seismogenic layer 
can be extracted from the data on QPOs for SGS. 
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Figure 1: Fractional frequency of nodeless spheroidal and torsional elastic oscillations as a 
function of fractional depth of seismogenic layer. 

2.2 Torsion mode 

For the torsional oscillations locked in the crust, the constants Ae and Be are eliminated 
from the following boundary conditions 

U^\r=R, = 0, U^\r=R = [(f>R X R]^, (17) 

0« = a.MV,P,(C). V,= (o.|.-I^|), (18) 

First is the no-slip condition on the core-crust interface, r = Rc, implying that the amplitude 
of differentially rotational oscillations is gradually decreasing from the surface to the core. 
The boundary condition on the star surface, r = i?, is dictated by symmetry of the general 
toroidal solution of the vector Laplace equation. The support of this last boundary condition 
lends further considerations showing that it leads to correct expression for the moment of 
inertia of a rigidly rotating star. The resultant algebraic equations steaming from above 
boundary conditions lead to 

TDe+2 

Ae = Are, B, = -U,Rl'^\ = i^^^+i - • (19) 



5 



Tedious calculation of integrals for Mt and Kt leads to 

47r£(£+l) pR^ 
(2£+ l)(2£ + 3) (1 - A2^+i)2 
(2£+l)2 
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(21) 
(22) 



In the limit of zero-size radius of the core, A = (Rc/R) 0, corresponding to torsional 
oscillations in the entire volume of the star we regain the early obtained spectral formula for 
the global nodeless torsional elastic vibrations I'iote) = [(2£ + 3)(2£ — 1)]"'^'^^ showing that 
in case of global torsional oscillations the lowest overtone is of quadrupole degree [SI [131 fTTj . 
However, this is not the case when we consider torsional nodeless oscillations locked in the 
seismogening layer of finite depth Ai? = R — Rc. For i = 1, equations (l20l) and (|2T|) are 
reduced to 



*V ' J 15(1 -A3)2 

\3 

Kt{i, A) = 871 fiR^ 



[1 - SA^ + 9A^ - SA^] 
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(23) 
(24) 
(25) 



In the limit when the core radius tends to zero, the stiffness Kt{i = 1, A = 0) ^ and the 
mass parameter getting the form of the moment of inertia of absolutely rigid solid star of 
mass A4 and radius R: Mt{i = 1, A = 0) = (2/5)A^i?2. This latter case corresponds to 
the rigid rotation. The above consideration again shows that the dipole overtone exhibits 
features of the Goldstone soft mode owing its emergence to the trapping of torsional shear 
oscillations in the peripheral crust of finite depth. 

The general spectral equation for the fractional frequency of nodeless torsional oscilla- 
tions of arbitrary multipole degree £, computed with aid of equations (|2UI) and ([21]), can be 
presented in the following analytic form 



[(£ + 2)(^-l)]p,(z/o,A) 
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(27) 



1 



(2£ 



A2^+i[(2£ - 1)(2£ + 3) - (2£ + 1)2A2 + (2^ + 3)A2^+i] 
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Figure 2: Theoretical curves for the frequency of spheroidal (dashed) and torsional 
(solid) nodeless elastic oscillations computed with aid of spectral formulae for frequency 
of spheroidal, eq.(15), and torsional, eqs. (26)-(27), modes as functions of multipole degree 
in juxtaposition with data (symbols) on QPOs for SGR 1900+14 and for SGR 1806-20. 

The usefulness of such representation is extensively discussed in|6j. The fractional frequency 
as a function oi h = AR/R is pictured in down panel of Fig.l which show that the lowest 
overtone of torsional vibrations trapped in the crust is of dipole degree and that dipole 
overtone of differentially rotational vibrations of the crust against core possesses properties 
of the Goldstone's soft mode. 

3 Application to SGR 1900+14 and SGR 1806-20 

The obtained spectral formulae (15) and (26)-(27) describe the frequencies of both spheroidal 
and torsional nodeless oscillations as functions of the multipole degree i. The natural unit of 
frequency z/q of shear elastic vibrations and the fractional depth h of peripheral seismogenic 
layer are input parameters carrying information about material properties of neutron star 
matter (density, shear modulus) and geometrical sizes of star and seismoactive zone. Con- 
sidering the observation data for SGR 1900+14 and SGR 1806-20, we demonstrate here how 
the obtained spectral equations can be used to eliminate some uncertainties in identification 
of QPOs. 

First, we examine the agreement of obtained spectral formula (l26!) - (l27j) for torsion mode 
with identification of the QPOs frequencies from interval 30 < z/ < 200 Hz with frequencies 
of nodeless torsional vibrations of multipole degree i from interval 2 < £ < 12 suggested 
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Figure 3: Theoretical predictions for the frequency of spheroidal (dashed) and torsional 
(solid) nodeless elastic oscillations as a function of multipole degree i in juxtaposition with 
data (symbols) on QPOs during the flare of SGR 1806-20. Based on the results of this 
latter work, our calculations suggest that low frequency QPOs, pictured by triangles, can be 
identified as dipole toroidal and dipole spheroidal nodeless vibration, respectively: i^(o^i) = 
18 Hz and u{oSi) = 26 Hz. 

in[3]. In so doing we use the proposed in these latter works identification of QPOs in SGR 
1900+14 data [namely, 1/(0^2) = 28 Hz; uioU) = 53; Hz //(o^e) = 84 Hz, i/(otii) = 155 Hz 
borrowed from Table 1 of paper[3]] as reference points and vary parameters z/q and h entering 
our spectral formula (15) for torsional mode so as to attain the best fit of these points. The 
result of this procedure is shown in upper panel of Fig. 2 by solid line. Then, making use of 
the fixed in the above manner parameters uq and h, we compute (with the aid of spectral 
formula ffT^ ) the frequency of spheroidal mode //(qS^). The application of this procedure to 
modal analysis of QPOs data for SGR 1806-20 is pictured in down panel of Fig. 2. Based 
on proposed in the above mentioned paper identification of the following points 1^(0^2) = 30 
Hz; i^(o^6) = 92 Hz and 2^(0^10) = 150 Hz we extract parameters z/g and h entering in our 
spectral formulae for torsional mode, equations (1261) - (1271) . 

In Fig. 3 set out very clearly the case for identifying two non-identified before low- 
frequency QPOs in data for SGR 1806-20 (Israel et al 2005), namely = 18 Hz and z/ = 26 
Hz, highlighted by triangles. 

4 Summary 

The obtained spectral formulae for the frequency of nodeless elastic vibrations trapped in the 
finite-depth seismogenic layer may be of some interest in its own right from the viewpoint of 
general theoretical seismology [9] in the sense that they can be utilized in the study of seismic 
vibrations of more wide class of solid celestial objects such as Earth-like planets. One of 
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the remarkable feature of considered model is that the dipole nodeless overtones possess 
properties of Goldstone soft modes, that is, the dipole overtones emerge if and on only if the 
elastic vibrations turn out to be locked in the peripheral layer of finite thickness. It is shown 
that obtained spectral equations are consistent with the existence treatment of low-frequency 
QPOs in the X-ray luminosity of flares SGR 1900+14 and SGR 1806-20 as caused by quake- 
induced torsional nodeless vibrations. What is newly disclosed here is that previously non- 
identified low-frequency QPOs in data for SGR 1806-20 can be attributed to nodeless dipole 
torsional and spheroidal vibrations, namely, i^(oti) = 18 Hz and t^(oSi) = 26 Hz. 
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